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ABSTRACT 

Context. The global size and spatial distribution of dust is an important ingredient in the structure and evolution of protoplanetary 
disks and in the formation of larger bodies, such as planetesimals. 

Aims. We aim to derive simple equations that explain the global evolution of the dust surface density profile and the upper limit of the 
grain size distribution and which can readily be used for further modeling or for interpreting of observational data. 
Methods. We have developed a simple model that follows the upper end of the dust size distribution and the evolution of the dust 
surface density profile. This model is calibrated with state-of-the-art simulations of dust evolution, which treat dust growth, fragmen- 
tation, and transport in viscously evolving gas disks. 

Results. We find very good agreement between the full dust-evolution code and the toy model presented in this paper. We derive 
analytical profiles that describe the dust-to-gas ratios and the dust surface density profiles well in protoplanetary disks, as well as the 
radial flux by solid material "rain out", which is crucial for triggering any gravity assisted formation of planetesimals.We show that 
fragmentation is the dominating effect in the inner regions of the disk leading to a dust surface density exponent of —1.5, while the 
outer regions at later times can become drift-dominated, yielding a dust surface density exponent of —0.75. Our results show that 
radial drift is not efficient in fragmenting dust grains. This supports the theory that small dust grains are resupplied by fragmentation 
due to the turbulent state of the disk. 

Key words, accretion, accretion disks - protoplanetary disks ~ stars: pre-main-sequence, circumstellar matter - planets and satellites: formation 



1. Introduction 

The dust content is fundamental for many aspects of the evo- 
lution, structure, and observations of protoplanetary disks. Not 
only is this the material out of which teiTestrial planets and 
the cores of giant planets form, but it also dominates the opac- 
ity, thereby determining the temperature structure and the ob- 
servational appearance of protoplanetary disks. The grain sur- 
faces also provide the ground for chemical surface reactions 
(e.g., Aikawa & Nomura 2006) and electron capture, so that dust 
strongly influences the ionization state and possibly the turbulent 
state of the disk (Sano et al. 2000; Bai & Goodman 2009; Turner 
et al. 2010). 

While there have been studies of gas line emission (e.g., 
Dutrey et al. 2007; Panic et al. 2008), most disk mass estimates 
are based on observations of dust emission, assuming a con- 
stant dust-to-gas ratio (see Andrews & Williams 2005, 2007). 
However, it is well known that the dust-to-gas ratio in disks 
should be an evolving quantity (e.g., Weidenschilling 1977a; 
Nakagawa et al. 1981; Kefler & Gail 2004; Brauer et al. 2007; 
Ciesla 2009). The problem with our understanding of dust trans- 
port is that it depends on the structure of the disk, which is un- 
certain, and also on the size distribution of the grains, which is 
also poorly understood. 

Several authors have investigated simplified models of dust 
transport and size evolution: for example, Garaud (2007), Kornet 
et al. (2001), Rozyczka et al. (2004), Klahr & Bodenheimer 
(2006), or Hughes & Armitage (2010). A more complete pic- 
ture is given by global simulations that self-consistently treat 



both the transport and collisional evolution, such as the works 
by Weidenschifling (1997), Brauer et al. (2008b), Bimstiel et al. 
(2010), or Okuzumi et al. (2011). The last models (which we 
will call full models as opposed to the simplified models, such 
as the one presented in this paper), however, are often so com- 
plicated that first it is difficult to derive a physical understanding 
from them, and second they are computationally expensive. 

We therefore put forward a model that is motivated by and 
based on the full simulations mentioned before (in this case, 
Birnstiel et al. 2010); however, it is very simplified so easier 
to understand, easier to implement, and computationally much 
less intensive. The aim of this work is to identify a few equa- 
tions that govern the spatial evolution of dust as found in the 
full simulations. We deiive upper limits for the size distiibution 
from which we can estimate the drift velocity of the largest par- 
ticles. The size distribution is then represented by just two grain 
sizes, one size representing the small dust that is coupled to the 
gas and another size representing the largest particles, which can 
be drifting inwards at a much higher velocity. While mathemati- 
cally similar, the basic approach of this work is thus quite dilfer- 
ent from other simplified models (e.g. Kornet et al. 2001 ; Garaud 
2007): we do not estimate the outcome of the physical processes, 
but instead identify the important concepts in the full models to 
construct a simplified model that is calibrated to reproduce the 
outcome of the full model. 

Klahr & Bodenheimer (2006) have estimated how much 
solid material can be concentrated in decimeter-sized radially 
drifting particles in local dust traps. Once these over-densities 
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Fig. 1: Dust surface density distributions as function of radius and grain size at 10 , 10 , and 3 x 10 years. In this simulation, the 
turbulence parameter was taken to be 10^^^. Overlayed are the representative sizes for a fragmentation limited size distribution 
(solid black lines) and for a drift limited size distribution (dashed black lines) according to Eq. 8 and Eq. 18, respectively. The 
dash-dotted line denotes the time-dependent estimate of the largest grain size a\(t) as defined by Eq. 28. 



reach Toomre unstable values and exceed the local Roche den- 
sity, planetesimal formation via self gravity in the dust layer is 
the consequence (Johansen et al. 2006). The latter mechanism 
has so far only be demonstrated for ad hoc local dust to gas ra- 
tios and size distributions. What is missing is a dust evolution 
model, as the one presented in this work, which predicts a local 
dust flux and the size distribution for models of gravoturbulent 
planetesimal formation (Klahr & Johansen 2008). 

The structure of this paper is as follows: in Section 2, we 
introduce a few concepts and formulae that are instrumental for 
the rest of this work. Section 3 presents some fiducial simulation 
results of the complete model that are used in Section 4 to derive 
and test a simple model of dust transport. Section 5 makes use of 
this model to discuss and explain some of the implication of it, 
such as the observational significance of the model, which pro- 
files of the dust surface density or the dust-to-gas ratio are to be 
expected under certain physical conditions, as well as the result- 
ing dust mass flux. Our findings are summarized in Section 6. 

2. Background 

Dust particles in protoplanetary disks are subject to radial 
drift (which depends on the pressure structure of the disk, see 
Weidenschilling 1977a), gas drag (i.e., particles following the 
gas accretion flow) and turbulent mixing. These effects depend 
on the size of the particles, which is also an evolving quan- 
tity because grains grow by sticking collisions. However there 
are several mechanisms that can limit collisional growth of dust 
particles: charging effects (Okuzumi 2009), bouncing (Guttler 
et al. 2010; Zsom et al. 2010), fragmentation and radial drift 
(Weidenschilling 1977a; Brauer et al. 2008a; Birnstiel et al. 

2010) . 

The charging and bouncing barriers are very much depen- 
dent on the material, porosity, or fractal dimension of the grains, 
and are still strongly debated topics (e.g., Wada et al. 2011). 
In one way or another these barriers need to be overcome in 
order to explain the existence of larger bodies such as aster- 
oids or planets. The fact that disks are observed to be rich in 
small dust particles (for a recent review, see Williams & Cieza 

2011) tells us that fragmentation must be effective in wide re- 
gions of the disk, otherwise coagulation would quickly render 
the disk optically thin even in the infrared wavelength regime 
(see Dullemond & Dominik 2004; Birnstiel et al. 2009). We take 
this observational fact as motivation for only considering frag- 



mentation/cratering and radial drift as growth-limiting mecha- 
nisms (Dominik & Dullemond 2008). For describing both of 
these effects, it is convenient to define the Stokes number St, 
which is the ratio of the stopping time of a particle fstop and the 
turn-over time of the largest eddy in the turbulent environment, 

St=^, (1) 

which under typical assumptions (Epstein drag law, compact 
spherical particles, isothermal gas density profile, and fL = ' , 
see Cuzzi et al. 2001) simplifies for particles near the disk mid- 
plane to 

St=^^. (2) 
2 

Here a is the radius and ps the internal density of the dust ag- 
gregate (taken to be 1.6 g cm^-') and Zg the gas surface density. 
The Stokes number is a dimensionless quantity that describes the 
aerodynamic properties of a dust particle and thus its coupling 
to the gas flow, i.e. particles with different sizes, densities, or 
shapes but the same Stokes number are aerodynamically identi- 
cal. 

Fragmentation of dust particles stops further growth because 
relative velocities tend to increase with the Stokes number of 
the particle (for St < 1). They can reach values well above the 
fragmentation threshold velocity Uf (Brauer et al. 2008a). Lab 
experiments measured threshold velocities of around 1 m s^' 
for the onset of fragmentation of silicate dust grains (e.g., Blum 
& Wurm 2008). It has been suspected that icy particles fragment 
only at higher velocities due to the ~10 times higher surface 
energies (Wada et al. 2008; Gundlach et al. 2011). Recent nu- 
merical studies of Wada et al. (2009) found that the fragmenta- 
tion velocity for <10;um sized icy aggregates could be as high as 
50 m s^'. However, experiments with silicate dust grains find a 
fragmentation velocity that decreases with grain size (Beitz et al. 
201 1), which can at least partly be attributed to the growing im- 
portance of inhomogeneities in larger grains (Geretshauser et al. 
201 1). In the following we discuss mostly results for a threshold 
velocity of 10 m s^', since we are interested in the global evolu- 
tion of dust, i.e. the regions beyond the ice line. The Stokes num- 
ber, at which the typical impact velocity of similar-sized grains 
equals the fragmentation threshold velocity can be approximated 
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Fig. 2: As in Fig. 1 but with a higher turbulence parameter at = 10 ^. 



by (see Birnstiel et al. 2009) 



Stf, 



1 Mf 



(3) 



where o-t is the turbulence parameter (Shakura & Sunyaev 1973) 
and Cs is the isothermal sound speed and we assume St < 1. 
The value of at is still largely unknown, but observations of ac- 
cretion signatures or disk lifetimes suggest values between lO^-' 
and 10^^. Theoretical models of the magnetorotational instabil- 
ity tend to reproduce these ranges of values (e.g. Johansen & 
Klahr 2005; Dzyurkevich et al. 2010), but global simulations 
also suggest a radial dependence of (see Sano et al. 2000; 
Flock et al. 2011). The simulations shown in this paper were 
carried out with radially constant values of at, but the model in- 
troduced in this work is able to treat also radial variations of a,. 

Apart from fragmentation, radial drift can also be a limiting 
factor for grain growth. Radial inward drift of dust particles oc- 
curs because the head wind induced by the sub-Keplerian gas 
disk removes angular momentum from the dust particles. The 
drift velocity (Weidenschilling 1977a) 



2 u„ 



St + St" 



(4) 



thus depends on the gas pressure gradient since Ujj is the dif- 
ference between the orbital velocity of the dust grain and 
the gas (i.e., the head wind velocity), which is given by (see 
Weidenschilling 1977a; Nakagawa et al. 1986) 



dP 1 



dr 2pgQk 



(5) 



Here P denotes the gas pressure, pg the mid-plane gas density, 
and the Keplerian angular velocity. It is straight-forward to 
also include reduction factors of the radial drift speed in Eq. 4, 
as found in Johansen et al. (2006), however this is not the subject 
of this work. 



3. Typical simulation results 

In this section we review some typical simulation results of the 
complete model that treats dust growth, fragmentation and cra- 
tering as well as radial transport due to gas drag, radial drift 
and turbulent mixing. These results are used in Section 4 to de- 
rive a simple model that reproduces the overall evolution of the 
dust surface density. The gas surface density is viscously evolv- 
ing using the typical turbulent viscosity prescription (Shakura & 



Sunyaev 1973). The temperature is estimated taking passive irra- 
diation and viscous heating into account. A detailed description 
of the model can be found in Birnstiel et al. (2010). Figures 1 
and 2 show snapshots from simulations of a 0. 1 Mq disk around 
a solar mass star. The contours denote cr(a, r), the dust surface 
density distribution per logarithmic size bin, which is defined 
such that 



f 

Ja=0 



cr(fl, r) dlna 



(6) 



is the total dust surface density. Throughout this paper we define 
r as the cylindrical distance to the central star. The initial gas 
surface density profile for the simulations is close to the self- 
similar solutions from Lynden-Bell & Pringle (1974) 



Sg(r)cx 



exp 



2-p 



(7) 



for a viscosity power-law index p of unity and a characteristic 
radius of 60 AU (see Fig. 3). The only significant deviation 
is in the inner parts of the disk where the viscous heating leads 
to a steeper temperature profile and consequently to a flatter sur- 
face density profile. Outside of about 3 AU, there is no deviation 
> 30% between the simulations and the self-similar solution. 
The evolution of the gas surface density for both simulations is 
shown in Fig. 3. 

Figures 1 and 2 differ only in the turbulence parameter at, 
which is 10~^ in the former and 10^^ in the latter case. There 
are a few fundamental features that are important for the further 
understanding of this work: 

- in the left panels of Fig. 1 and 2 it can be seen that grains 
in the inner regions quickly grow and reach a maximum 
size, the so-called fragmentation barrier, which is discussed 
in Section 4.1.1. The size limit set by fragmentation is shown 
as black solid lines' . A steady state is reached in which frag- 
mentation and coagulation balance each other. We call these 
regions fragmentation limited. As we go to larger radii in the 
first panel, we see that grains at about 20 AU and beyond are 
still at the state of initial growth, as discussed in Section 4.3. 
In the central panel of Fig. 2 we see a similar behavior, but 
growth does happen on slightly shorter time scales because 
the growth time scale of small particles (St < at) is shorter 
due to the higher turbulent velocities. It is only for larger 



' More precisely, the solid line is the representative size in a fragmen- 
tation limited distribution that is by definition slightly below the upper 
end of the size distribution, as can be seen in Fig. 1. 
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particles that the growth time scale does not depend on the 
turbulence strength (Brauer et al. 2008a). 

- In the central panels of Figs. 1 and 2, we can see that grains 
also in the outer regions have grown significantly. The size 
distribution in the inner regions are basically unchanged 
compared to earlier times since they are still in coagulation- 
fragmentation equilibrium. For reasons which are discussed 
later, the barriers to growth, shown as dashed and solid black 
lines in Figs. 1 and 2 also evolve with time. Most signifi- 
cant is the fact that the drift limit (dashed lines, discussed in 
Section 4.1.2) is moving down to smaller sizes. 

- In the outer regions at later times, the drift barrier (dashed 
line) has moved below the fragmentation barrier (solid line) 
as can be seen in the central and right panels of Fig. 1. 
This means that dust particles in those regions are no longer 
limited by fragmentation, but instead they are more rapidly 
drifting away than being replenished by growth from smaller 
sizes. We call those regions drift limited. 

- Due to the higher strength of turbulence in Fig. 2, the frag- 
mentation limit shifted to smaller sizes compared to Fig. 1, 
which renders the disk fragmentation limited up to about 
70 AU, while in Fig. 1, only the inner regions up to about 
2 AU are fragmentation-dominated. 



4. The two-population model 

In this section, we develop an analytical understanding of dust 
evolution in size and spatial dimension. The simulation results 
of the comprehensive model, which have been introduced in the 
previous section are used to derive simple equations that char- 
acterize the size distribution. This enables us to develop a sim- 
ple model that can reproduce the evolution of the dust surface 
density of the full-scale simulations well. While the following 
section comprehensively explains the derivation of the model. 
Appendix B summarizes the model in a concise algorithmic 
form. 



4.1. Size limits 

We crudely simplify the grain size distribution by only two rep- 
resentative sizes, fl()(r) and fli(r), where ao{r) is the monomer 
size and oi (r) is the representative size of the large grain popu- 
lation (see also the sketch in Fig 4). While ao{r) is taken to be 
constant in time and space, oi (r) increases with time as particles 
grow and the disk structure evolves. This growth is not without 
limits. As explained in Section 2, we consider the two most im- 
portant barriers towards further growth, namely grain fragmen- 
tation and radial drift. We find that oi (r) can also decrease with 
time if the growth limits drop to smaller sizes. 



4.1.1. Turbulent fragmentation 

In the case where the fragmentation of dust grains due to turbu- 
lent relative motion is the limiting effect, analytical models of 
the grain size distribution have been derived in Birnstiel et al. 
(2011). In these stationary and local grain size distributions, 
a significant fraction of the dust mass is concentrated in the 
largest sizes, slightly below the limiting Stokes number Stfrac,. We 
parametrize this offset by an order-of-unity constant /f such that 
the representative size for the largest grains in a fragmentation- 
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Fig. 3: Evolution of the gas surface density profiles of the low 
turbulent (top panel) and the high turbulent simulation (bottom 
panel). The solid, dashed, and dash-dotted lines correspond to 
simulation times of 10"*, 10^, and 3 x 10*', respectively. 
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Fig. 4: Sketch of the two-population representation of the size 
distribution. The solid black curve is an arbitrary dust size dis- 
tribution with the upper limit of Ofrag. The dark shaded area is the 
small population, i.e. the part of the size distribution which is not 
influenced by drift velocities because particles of these sizes are 
low enough to be tightly coupled to the gas. The light gray area is 
the large population, i.e. the grains that are significantly drifting 
inwards. The total mass contained in the small and large popu- 
lation is given by Eq and Ei, while the representative size is ao 
and «!, respectively. 



dominated distribution can be stated as (using Eqns. 1 and 3) 

2 Sa M? 



/f 



3 n p^ai cl ' 



(8) 
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4.1.2. Radial drift 



4.1.3. Fragmentation by relative drift velocities 



Radial drift can also set a limit to the local size of dust grains if 
the time within which particles are removed by drift is similar 
or less than the time scale at which the particles are replenished, 
i.e. the growth time scale. To derive the growth time scale 

a 

^tjrow ~ 5 (9) 

a 

we follow Brauer et al. (2008a) by using the growth rate of 
monodisperse coagulation 

^ = (10) 
df ps 

and the approximate relative velocities between the similar-sized 
dust grains for turbulent motion (see Ormel & Cuzzi 2007) 



Am 



'3 fft StCs- 



(11) 



Assuming St < 1, the Epstein drag law, and a scale height of the 
dust distribution of (Youdin & Lithwick 2007) 



Ha = Ho 



St' 



(12) 



(13) 



the growth time scale can be written as 

1 

Tgrow — 7^r~^ 

where e = Sj/Sg is the vertically integrated dust-to-gas mass 
ratio and we neglected a factor of order unity. 

It should be noted that in this case (i.e. turbulent relative ve- 
locities), the growth time scale is independent of the particle size 
because the rate of change of the particle radius becomes directly 
proportional to the particle radius itself. The drift time scale 

Tdrift = — (14) 

depends on the drift velocity (cf. Eq. 4) and thus the Stokes num- 
ber of the particle. It is given by 



Tdrift 



St cf 



r 



where is the Keplerian velocity and 
^ dlnr 



(15) 



(16) 



is the absolute value of the power-law index of the gas pressure 
profile. By equating the growth time scale and the size dependent 
drift time scale, we derive an estimate of the maximum Stokes 
number that can be reached before radial drift removes the dust 
particles, 



Std, 



rift 



y2 



(17) 



or in terms of particle size (again assuming Epstein drag law), 

2 2d -1 ox 
Adrift = /d J 7 ■ (18) 

This boundary is not as sharp as the fragmentation barrier, as 
can be seen in Fig. 1. We have introduced another order-of -unity 
factor /d, which we calibrate later with the detailed numerical 
simulations. 

Strictly speaking, Eqns. 11 and 12 are not valid for very 
small particles that are very tightly coupled to the gas. However 
this limit was derived for drifting particles, which are per defini- 
tion only marginally coupled to the gas thus Eqns. 1 1 and 12 are 
valid approximations. 



So far we have only considered fragmentation by turbulent rela- 
tive velocities. However in the case of very low turbulent veloc- 
ities, one might expect that the relative drift speed of the parti- 
cles, either in the vertical or in the radial dimension, is still high 
enough to cause fragmentation. Dust grains in the tenuous disk 
atmosphere settle towards the mid-plane with velocities that can 
exceed the radial ones (Nakagawa et al. 1986). In a turbulent 
disk, mixing counteracts this downward motion and an equilib- 
rium state is reached on time scales less than 10^ years, which 
means that the vertical motion is not significant throughout the 
lifetime of the disk. In contrast to this, radial drift motion does 
not vanish towards the mid-plane and is therefore a systematic 
velocity that stays relevant throughout the lifetime of the disk. 
We therefore only focus on the fragmenting effects of the radial 
dust motion. We show in Section 5.1 that this effect is irrelevant 
for the simulations of a turbulent disk. For laminar disks, this 
effect can become relevant. 

To derive a size limit for this case, we need an estimate of 
the relative velocities induced by radial drift. From Eq. 4, we can 
derive the relative drift speed between two particles for Stokes 
numbers below unity 



Am, 



drift 



r |sti -st2| 



(19) 



Expressing the Stokes number of the smaller particle as St2 = 
• Sti, and equating Eq. 19 with the fragmentation threshold 
velocity Mf, we derive 



Stdf 



Mf Vk 



ycl{\-Ny 



(20) 



which is the largest Stokes number that can be reached before 
the relative drift speed exceeds the fragmentation threshold. The 
most important collision types determining the mass gain/loss of 
the largest bodies is typically collisions with similar sized bod- 
ies. A^ = 0.5 is therefore a reasonable assumption. 

4.2. Mass fractions & transport equation 

In order to represent the size distribution of the dust with only 
two representative sizes, we also need to know the ratio of mass 
contained in large and small particles. We represent the mass 
ratio between the two populations by another yet unknown factor 

fm-> 



2i(r) = EdW-A,(r) 
2o(r) = IdM-(l-/„(r)), 



(21) 
(22) 



where is the surface density in the small grains while is the 
surface density in large grains. This representation is schemati- 
cally shown in Fig. 4. In the case of a dust size distribution which 
is in coagulation/fragmentation equilibrium, small dust is replen- 
ished. In the case of a drift-limited size distribution, fragmenta- 
tion is not efficient enough since particles drift inwards before 
reaching sizes at which they fragment. 

For particles with a Stokes number below unity, this allows 
us to evolve the dust distribution by a single advection-diffusion 
equation 



5Sd 1 d 
dt r dr 



d 

r ["Lau — Dgas Eg — 



0, (23) 
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- Table 1 : Best-fit model parameters 
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Fig. 5: Grain size distributions (top panel) and the cumulative 
flux (bottom panel) for the simulation result which is displayed 
in the right panel of Fig. 1 at 1, 10, and 100 AU. The differ- 
ent radii are represented by different gray scale colors and the 
dashed lines mark the representative grain size oi . 




10-' 10- 

Stokes number 



Fig. 6: The dust flux distribution cr{a) ■u{a) as function of Stokes 
number for the same dust distributions as shown in Fig. 5. 



parameter 


drift limited 


fragmentation 




case 


limited case 


/f 




0.37 


/d 


0.55 






0.97 


0.75 



Notes. The values in this table represent the values of the model param- 
eters that best reproduce the results of the full simulations. 



4.3. Particle growth 

The mass flux and thus the whole dust evolution is strongly in- 
fluenced by the growth time scale. This is due to the fact that 
the dust-to-gas coupling and thus the radial velocity strongly 
depends on the particle size. Thus, particles in the outer parts, 
which grow less quickly, start to drift at later times, once they 
have reached sizes large enough to be influenced by drift. From 
the growth time scale (i.e. the size-doubling time) given in Eq. 13 
we can estimate the time it takes to grow from the smallest size 
ao to a maximum size ai 



^grow 



'^grow * 111 



where 



fll = mm (^Cldrift, flfrag 



(26) 



(27) 



In order to take this delayed drift into account, we define a time 
dependence of the larger representative size 



fli (f) = min 



ai, flo • exp 



(28) 



using the mass weighted velocity 



This is an oversimplified estimate since different growth mech- 
anisms are at play at different sizes, however we found a good 
agreement between our estimate and the simulation result. This 
can be seen in Figs. 1 and 2, where the dash-dotted curve rep- 
resents the time dependent upper size limit according to Eq. 28. 
Only at early times (left panels of Figs. 1 and 2), fli(f) differs 
from the fragmentation or diift baiTier 

It can also be seen that this estimate is slightly too low, i.e. it 
is slightly below the upper size of the simulated dust distiibution 
at all radii > 20 AU in the left panel of Fig. 2. This is due to 
the fact that the growth time scale of Eq. 13 implicitly assumes 
growth to be driven only by equal sized collisions from turbu- 
lent motion according to Eq. 1 1 . However, the initial stages of 
particle growth involve different turbulent regimes (see Ormel & 
Cuzzi 2007), which lead to faster growth for higher values of at. 



M= (l-/„(r))- MO -Ml, (24) 

and the gas diffusivity Dga,. Here, mq and mi denote the radial 
velocity of the dust sizes aq and oi, respectively. The radial ve- 
locity of dust grains is the sum of the drift velocity given by Eq. 4 
and the gas drag induced velocity 

1 

i ~\~ ot 

where Mgas is the gas radial velocity. A more thorough justifica- 
tion of Eq. 23 is given in Appendix A. 



4.4. Calibration and test cases 

In the semi-analytic model described above, we have introduced 
three factors /f , /d, and /„. These factors are necessary since our 
model is not derived in a rigorously analytical way, but includes 
several assumptions and estimates. In order to still achieve a 
good agreement with the detailed numerical simulations, we 
need to calibrate these order-of-unity factors by comparison to 
a grid of the full simulation results. 

In order to do this, we divide the grain size distribution at 
each radius in two populations: the large population, which are 
all grains whose drift velocity is higher than the gas drag in- 
duced velocity, while everything that is smaller belongs to the 
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at 5.4 AU 



small population. The gas radial velocity can be written as (see 
Takeuchi & Lin 2002) 



3 fftCj 



{2-p~q), 



(29) 



where p is defined as l.gCCr^'' and q as Tccr^'', and is taken 
to be constant. Thus, the gas drag velocity is lower than the drift 
velocity for particles with a Stokes number larger than 



Ste 



3\2-p-q\ 
> — at. 



Irl 



0.5, r 



(30) 



which for typical values (p = 1, q 
yields Stgq > Q't/2. 

As described in Section 4.1, the fragmentation is a rather 
strict upper Umit to the size distribution. The dust mass flow is 
therefore not dominated by grain sizes at the fragmentation bar- 
rier but slightly below it. We therefore calculated a flux-averaged 
grain size for the large population 



a{r 



where 



(r{a, r) u{a) da 




cr{a, r) u{a) 



da. 



St 



•eq 



22g 



(31) 



(32) 



is the grain size corresponding to a Stokes number of Stgq and 
u{a) is the total radial velocity of grain size a. We compared the 
outcome of the simplified model to a grid of 39 full simulations." 
Good agreement of this flux averaged size with Ofrag was found 
by tuning /f to a value of 0.37. 

For the case of a drift limited size distribution, the size limit 
is not as strict as in the fragmentation limited case, i.e. there can 
exist grains slightly larger than adrift- The same method as above 
yields therefore a value of /d of 0.55, which is higher than in the 
case of a fragmentation limited size distribution. These values 
of /f and /d should thus be universal, as long as the collision 
outcome model is not changed. 

The dashed and solid lines in Fig. 1 denote the represen- 
tative grain size for a drift limited and a fragmentation limited 
maximum grain size, respectively. It can be seen that these sizes, 
together with the time dependence introduced in Section 4.3, de- 
scribe the limits of the grain size distribution well. 

In order to properly capture the radial transport of dust with 
this simple model, we need also to calibrate the parameter /„, 
which describes how the mass is distributed between the two 
populations. Typical size distributions can be seen in the top 
panel of Fig. 5, denoted by the solid lines. Both the drift and 
the fragmentation limited distributions have most of the mass in 
the large grain population. However fragmentation limited dis- 
tributions (e.g., the solid black line in the top panel of Fig. 5) 
have a larger fraction of mass in the small sizes compared to 
drift limited distributions (solid gray lines). The value of /,„ thus 
depends on whether drift or fragmentation is the limiting factor 
We found typical values of /m around 0.75 for the fragmentation 
limited case and around 0.97 for the drift limited case. All values 
are summarized in Table 1 . 



^ We carried out 39 simulations using as initial condition either 
self-similar solutions (Eq. 7) with = 20, 60, 200 AU, at = 
10^^, 10^^, 10^^, disk masses of 0.01 and 0.1 Mq, fragmentation ve- 
locities of 1, 3, and 10 m s^' or power-law models with gas surface 
density exponents of -1 with outer radii of 250 AU. 
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at = 10- 
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- two-populations model 

gas flux X dust-to-gas ratio 
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Fig. 7: Total inward mass flux of dust for the full simulations 
(solid line), for the two-population model (dashed line) and the 
scaled gas accretion rate (dot-dashed line) in earth masses per 
year at 5.4 AU. The upper and lower panel correspond to the 
weak and strong turbulence simulations, respectively. 



It should be noted that the flux is typically determined by 
the large population. This is shown in Fig. 6. The plotted quan- 
tity cr(a) • u(a) is defined such that integration over log(a) 
(or equally log (St)) yields the total dust flux. This means that 
a one-population model, which neglects the small population 
can describe the dust evolution in the drift limited case well. 
However in the case of a fragmentation-dominated distribution 
or an outward-spreading gas disk, the small population can con- 
tain enough mass or have a significant velocity to be relevant and 
generally better agreement is obtained if this population is taken 
into account, while the additional complications and computa- 
tional costs of the two-population model are minimal. 

Another important test case is whether the total dust mass 
flux is reproduced. In Fig. 7 we compare the total dust flux from 
the full simulations (solid line) to our two-population model 
(dashed line). The total dust flux is calculated as 



r 

Jo 



Md(r) = 2;7rr cr(a, r) M(a) dlna. 



(33) 



The dash-dotted Une in Fig. 7 represents the gas accretion rate 
multiplied by the initial dust-to-gas ratio, i.e. the dust mass flux 
that would be found if all dust particles were perfectly coupled 
to the gas accretion flow. The error of the toy model is always 
within a factor of two, which is a surprisingly good result con- 
sidering the simplicity of the model. The results are roughly 
consistent with the order-of-magnitude estimates of Klahr & 
Bodenheimer (2006). 

Finally, Fig. 8 compares the profiles of the dust-to-gas ratio. 
It can be seen that the two-population model reproduces the full 
simulation results reasonabl well. 
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Fig. 8: Evolution of the vertically integrated dust-to-gas ratio for the weak and the strong turbulence simulation (upper and lower 
panels respectively). The panels on the left hand side show the results of the full simulation, the panels on the right hand side those of 
the simplified model. The dotted line is the initial condition which is taken to be a constant dust-to-gas ratio of 0.01 up to 180 AU. 
Even the smallest grains at the low gas densities beyond 180 AU would decouple and quickly drift inside. We chose to remove 
the dust beyond 180 AU to avoid a wave of dust sweeping inwards, which would be caused by an initial condition with constant 
dust-to-gas ratio throughout the disk. 



5. Discussion 



Stdrift and Stdf , which is given by 



The two-population model presented in Section 4 was shown 
to reproduce the time evolution of several important quantities, 
such as the dust-to-gas ratio, the dust mass flux and also the size 
of the drifting particles. In this section we discuss some of the 
applications of this model. 

5.1. Fragmentation by drift velocities 

In the following subsection we investigate whether the relative 
velocities due to radial drift of the particles can become high 
enough to play a role. For this to be the case, the size limit of 
Eq. 20 must be smaller than the limits set by both turbulent frag- 
mentation (cf., Eq. 3) and by transport (cf., Eq. 17). These two 
condition can be rewritten as 



Vk 3 



and 



£ > 



Uf 



(34) 



(35) 



respectively. Considering a typical disk where y ~ 2.75 and as- 
suming = 0.5, we find that for strong turbulence (at > 10^^), 
the first condition is never fulfilled even for high fragmenta- 
tion velocities of ~ 10 m s^'. For weaker turbulence such as 
at < 10^^, the regions beyond ~ 10 AU can fulfill Eq. 34. 
Eq. 35 sets constraints from the opposite direction, such that the 
region between about 10 and 60 AU could be influenced by drift 
induced fragmentation. 

To quantify the error we do by considering only the size lim- 
its by drift and fragmentation, we consider the ratio between 



Stdri 



Stdf 



2.5 • 



0.01 



Uf 



10ms- 



-r' f^r^ 

ly vioAu/ 



(36) 



It is important to note that this deviation only needs to be consid- 
ered in regions where it is larger than unity and turbulent frag- 
mentation does not play a role. Since the regions further inside 
are dominated by turbulent fragmentation, the largest error we 
do is overestimating the upper end of the distribution by a fac- 
tor of about 2.5. As we go to 60 AU the error is decreasing while 
the regions beyond 60 AU are in any case limited by particle drift 
and not by fragmentation. Decreasing the fragmentation thresh- 
old velocity seems to increase the error, however in this case 
the region where drift induced fragmentation applies, then dis- 
appears because turbulent fragmentation becomes more impor- 
tant. 

Most importantly, we see that Eq. 36 depends on the dust-to- 
gas ratio. Since the region beyond 10 AU is strongly drifting, the 
dust-to-gas ratio e is quickly decreasing as can be seen in Fig. 8. 
Therefore, we can safely ignore the size limit set by Eq. 20 if 
we are only concerned about the evolution of the dust surface 
density or the upper end of the dust size distribution. 

We also carried out simulations with a very small amount 
of turbulence (a^ = 10~^). As expected, particles grow to large 
Stokes numbers until fragmentation due to radial drift sets in. 
However even in this case, the dust-to-gas ratio drops so quickly 
that after 1 Myr the drift limit (Eq. 17) becomes everywhere 
smaller than the limit set by drift induced fragmentation (Eq. 20), 
such that grains do not fragment anymore. In terms of the pres- 
ence of small dust, this effect does matter: it means that a drift- 
limited size distribution cannot sustain efficient fragmentation 
on long timescales and thus contains much smaller amounts of 
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small dust grains (see right panel of Fig. 1). The observational 
fact that disks are observed to be rich in small dust for several 
millions of years thus seems to be in contradiction with a com- 
pletely drift-dominated disk unless there is an external source of 
small dust (see Dominik & Dullemond 2008) or small dust is 
mixed and fragmented into surface layers with stronger turbu- 
lence, an effect that we cannot take into account in a vertically 
integrated dust evolution code. We have shown that radial drift 
alone cannot sustain the relative velocities to keep fragmenta- 
tion active for the observed lifetimes of protoplanetary disks. 
Fragmentation together with vertical mixing is therefore likely 
driven by turbulent motion. 

5.2. Analytical dust surface density profiles 

The evolution of the dust-to-gas ratio is important for theories 
of planetesimal formation (e.g., Johansen et al. 2007, 2011) and 
planet formation (see Klahr & Bodenheimer 2006; Lissauer & 
Stevenson 2007; Mordasini et al. 2009). Also a constant dust- 
to-gas ratio is usually assumed for deriving disk masses from 
dust continuum observations (e.g., Andrews & Williams 2005, 
2007), which is one of the biggest sources of error in disk mass 
estimates. In this section, we investigate how the dust surface 
density and thus the dust-to-gas ratio changes with time for a 
typical model of a circumstellar disk, starting with a constant 
dust-to-gas ratio of 0.01. The following discussion is based on 
two important aspects: 

- Most of the (dust) mass resides in the outer regions of the 
disk. Due to the long evolution time scale at these radii, these 
regions provide a mass reservoir for the inner parts of the 
disk, as was also found in Kornet et al. (2004) or Garaud 
(2007). 

- As discussed in previous sections, the dust mass is concen- 
trated towards the upper end of the size spectrum, which is 
also where the drift velocity is the highest. This means that 
the dust flux is governed by the upper end of the size distri- 
bution. 

A semi-steady dust surface density is then set by the rate at 
which the outer regions provide dust which is inward drifting 
with a velocity of u. Mass conservation then dictates that the 
dust accretion rate 

Md = 27rrEdM (37) 

has to be constant for all r. This yields a dust surface density 
profile 

^ Md 1 
In ru 

For a drift-dominated distribution, we can derive the drift veloc- 
ity of the representative size fli(f) from Eqns. 17 and 4, which 
results in a surface density profile given by 



(38) 



-d.drift 




(39) 



of St < at/2, see Eqns. 30 and 32) are coupled to the gas, while 
for somewhat larger particles, radial drift starts to play a role. For 
a dust surface density distribution with an upper size limit larger 
than fleq (cf. Eq. 32), we can derive a stationary surface density 
profile as we did before but using Eq. 8 as upper size limit. This 
yields 



■^d,frag 



3Md Ct 



cx- 



271 ffJU^ 



(40) 



which for constant values of at and Uf is proportional to r 
as the Minimum Mass Solar Nebula profile (Weidenschilling 
1977b; Hayashi 1981). 
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Fig. 9: Comparison between simulation results and the simple 
estimates derived in Section 5.2 after 3 Myrs of dust evolution. 
The upper panel is for a turbulence value of ai = 10^^, the 
lower one for = 10^^. The estimates have been calculated 
using the dust accretion rate of the simulations at 50 AU which 
are 1.6 x 10^^ and 2.4 x 10^^M©/yr, respectively. The division 
between the different fit formulas is described in Section 5.2. 



A dust distribution that consists entirely of small dust grains 
(St < at/2 for the largest grains) should be so well coupled to 
the gas that the dust-to-gas ratio is constant, i.e. 



■-d,mixed 



Md Q. 



1 



27T3atcl \2-p-q\ 



(41) 



which for a gas surface density profile with index /? = 1 is pro- 
portional to r~^^^. 

For the case of a fragmentation-dominated distribution, such 
a result cannot be obtained as uniquely, because very small par- 



because the dust is just co-moving with the gas flow. 

We have now derived three different shapes of the dust sur- 
face density distribution. In Fig. 9, we compare our simple es- 
timates to the previously discussed simulation results of the full 



tides which are smaller than Oeq (i-e., they have a Stokes number model. Both simulations are for a 0.1 Mq disk with = 60 AU 
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around a solar mass star assuming a fragmentation velocity of 
10 m s^'. The turbulence in the simulation of the upper panel 
is weaker (o-t = 10~^) than in the bottom panel (ct = 10~^). 
It can be seen that the simulation in the upper panel becomes 
drift-dominated in the outer parts (i.e. ajrift < flfrag)- The drift 
estimate (cf., Eq. 39) agrees well with the simulation result in 
the regions between about 2 and 200 AU, while the inner parts 
follow the fragmentation limited estimate, Eq. 40. The estimates 
do not capture the drop-off in the outermost parts of the disks 
because they assume a constant accretion rate, while in reality, 
the surface density and thus the accretion rate must decline to- 
wards the outer end of the disk. The stronger turbulence in the 
bottom panel causes the upper limit of the dust size distribution 
to decrease. The fragmentation-dominated region now reaches 
out to larger radii at all times while grains in the innermost re- 
gions become so small that they are well coupled to the gas and 
thus can be described by the mixing estimate of Eq. 41 . Only the 
regions beyond about 60 AU stay drift-dominated. For a lower 
fragmentation threshold velocity of 1 m s^', the entire disk was 
found to be fragmentation or mixing dominated. 

The results shown in Fig. 9 are for 3 Myrs of dust evolution. 
At earlier times, for example 1 Myr (not shown), the dust-to-gas 
ratio in the outer parts is still higher and thus the drift limited 
region is smaller, covering only the region beyond 4 and 150 AU 
for the low and high turbulence case, respectively. 

Andrews et al. (2012) have recently found a discrepancy in 
the sizes and shapes of the dust and gas surface densities in the 
disk of TW Hya, based on observations of 870 fim dust emis- 
sion and CO J 3-2 emission. In Fig. 10, we plot these observa- 
tionally derived profiles along with the analytical, drift limited 
dust surface density (using the observational gas surface den- 
sity in Eq. 39). The analytical profile has been scaled to fit the 
value of the dust profile from Andrews et al. (2012) at 10 AU. 
Such a profile would be expected for a moderate level of turbu- 
lence at late stages of disk evolution (the age of TW Hya is esti- 
mated to be about 10 Myrs, see Kastner et al. 1997; Webb et al. 
1999). Fig. 10 shows that the shape of the analytical solution 
agrees very well with the observations. One notable difference is 
the shape of the outer edge. Andrews et al. (2012) note that the 
remnant mismatch between their best-fit model and the 870 /im 
emission is likely related to the shape of the outer edge, which 
in their case is a sharp cut-off. In contrast to that, our model as- 
sumes a constant dust mass accretion rate through the disk. As 
stated earlier, the mass accretion rate has to drop off at some 
radius as the dust reservoir in the outer regions gets depleted. 
The two dust profile in Fig. 10 thus most probably represent two 
extreme cases, i.e. a too sharp and a too smooth outer edge. 

5.3. Evolution of the dust-to-gas ratio and dust accretion rate 

The results presented in the previous section can directly be 
translated into a profile of the dust-to-gas ratio. The absolute 
value of this profile however depends on the remnant dust ac- 
cretion rate, which is initially higher due to radial drift, but 
more quickly decreasing than the gas accretion rate, as shown 
in Fig. 7. The fact that the dust accretion rate is higher than the 
gas accretion rate multiplied by the initial dust-to-gas ratio is im- 
portant for theories of planetesimal or planetary core formation 
because this allows a dust trap such as a vortex or zonal flow to 
achieve local dust over-densities. Depending on the trapping effi- 
ciency and the rate at which dust is accumulating in the trap, the 
dust over-densities can reach critical values to trigger the stream- 
ing instability and gravoturbulent formation of planetesimals, as 




10° 10^ 

r [AU] 



Fig. 10: The best-fit models of Andrews et al. (2012) for the 
dust and gas profiles of TW Hya (solid and dashed black lines, 
respectively) and the analytical, drift limited dust profile as in 
Eq. 39 (solid grey line), using the gas profile of Andrews et al. 
(2012) is used as input. 



discussed in Klahr & Bodenheimer (2006) and Johansen et al. 
(2006). 

In the low turbulence case (top panels of Fig. 8) it can be 
seen that the disk regions inward of 4 AU get strongly enriched 
in dust because grains quickly grow to larger sizes and drift. This 
leads to a pile-up effect similar to Youdin & Chiang (2004). This 
behavior can be understood by dividing Eq. 40 by the gas sur- 
face density, which results in a dust-to-gas ratio proportional to 
5jjj(;g steady state a disks obey p + q = 3/2, we can also 

write 

efragCCT. (42) 

Thus the steep increase in temperature due to viscous heating 
causes the increase in the dust-to-gas ratio through the tempera- 
ture dependence of the gas disk profile. 

In the high turbulent case (bottom panels of Fig. 8), the dust 
enhancement is much less effective because of the more effec- 
tive fragmentation of dust grains due to higher turbulent veloci- 
ties. This causes the drift velocities to be lower than in the low- 
turbulent case and for a given accretion rate, this yields only a 
modest enhancement of the dust-to-gas ratio. 

It is important to note that so far we have considered the ver- 
tically integrated dust-to-gas ratio, which is what future obser- 
vations may be able to constrain. For theoretical works, the mid- 
plane dust-to-gas ratio is usually more important. Due to set- 
tling, this will be even larger than the integrated value, if grains 
are large enough to efficiently settle down to the mid-plane. For 
particles with Stokes numbers at < St < 1, Eq. 12 describes 
the dust scale height as function of Stokes number reasonably 
well. For the low turbulence case, the particles at the drift limit 
(Eq. 17) reach Stokes number of up to 0.1 at early times, later 
dropping < 0.01 (see also Fig. 5). According to Eq. 12, this 
means that the dust scale height and thus the mid-plane dust- 
to-gas ratio is enhanced by a factor of 3 to 10 compared to the 
vertically integrated values shown in Fig. 8. We find that at early 
times - up to 3 x 10^ years - the mid-plane dust-to-gas ratio 
in the whole disk is increased at least by a factor of 2 over the 
canonical values. After that time, the loss of dust mass coun- 
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teracts the settling effects. This result obviously depends on the 
initial condition of the dust-to-gas ratio, which in our case was 
taken to be a constant value of 0.01. 

6. Summary and conclusions 

In this paper, we have derived a simple model that describes the 
radial evolution of the dust surface density under the combined 
influence of growth and fragmentation of compact grains as well 
as radial transport mechanisms. This has been achieved by find- 
ing upper limits to the grain size distribution, which are func- 
tions of time and the physical conditions at each radius. Good 
agreement between this very simple numerical model and a 
full-fledged, computationally intensive dust evolution code was 
found. 

Important parameters are the fragmentation threshold veloc- 
ity Mf, the level of turbulence o-t, the initial dust-to-gas ratio, and 
the temperature and density profile of the gas disk. The effective 
dust transport velocity can be estimated by representing the dust 
distribution by only two characteristic grain sizes, the small and 
the large population. This allows us to derive analytical solutions 
for the dust surface density profiles in protoplanetary disks. 

Our findings are summarized in the following: 

- The simple two-population model presented in this work de- 
scribes the evolution of the dust-surface density and the evo- 
lution of the largest grain sizes well. Good agreement be- 
tween this simplified model and a full-fledged dust evolution 
code was found. Additionally, the dust-to-gas ratio, the dust 
flux, and the size of the drifting particles can be derived from 
it. 

- The upper end of the grain size distribution can be described 
as limited by turbulent fragmentation (cf., Eq. 8) or by ra- 
dial drift (cf., Eq. 18). Fragmentation due to relative radial 
drift alone (Eq. 20) is found to be ineffective in non-laminar 
disks (see Section 5.1). This supports the theory that disks 
are turbulent because radial drift alone cannot cause efficient 
fragmentation over long enough timescales to agree with ob- 
servations. 

- We derived three diff'erent analytical profiles of the dust- 
surface density for different physical conditions: firstly, for 
a strongly drifting dust distribution (Eq. 39), secondly for a 
fragmentation limited distribution (40), and thirdly for a dis- 
tribution where all grains are so small that they are well cou- 
pled to the gas (Eq. 41). The free parameter of the profiles 
is the dust accretion rate provided by the outer regions. The 
analytical profiles were found to fit to the simulation results 
of the full dust evolution code of Birnstiel et al. (2010) very 
well. 

- We found that at late times of their evolution, disks can be- 
come drift-dominated, which for typical gas disk profiles 
(SgXr^', Tcxr""^) leads to a dust surface density profile 
proportional to r""-^^. These results agree with the best-fit 
models for dust and gas profiles of the ~10 Myr old TW 
Hya disk, as recently found by Andrews et al. (2012). 

- In the case of a fragmentation limited distribution with rel- 
atively large grains present, the dust surface density profile 
becomes proportional to r^'^ as in the Minimum Mass Solar 
Nebula (see Weidenschifling 1977b; Hayashi 1981). 

- Similar to Youdin & Chiang (2004), the vertically integrated 
dust-to-gas ratio is strongly enhanced in the innermost re- 
gions if the dust is significantly drifting in the outer region 
of the disk. This pile-up is aided by the combined effects of 
drift and fragmentation. 



- An enhancement of the mid-plane dust-to-gas ratio was 
found in the low-turbulent simulation (o-t = 10~^). This is 
because particles can reach large enough Stokes numbers (up 
to 0.1) to efficiently settle to the mid-plane. In the case of 
efficient fragmentation, particles remain too small to signifi- 
cantly enhance the mid-plane dust-to-gas ratio. 

- As an important parameter for models of planetesimal or 
planetary core formation the radial mass flux in solids has 
been determined as a time and space dependent function. 
Depending on the disk parameters, the dust accretion rates 
can lie anywhere from a factor of a few up to two orders 
of magnitude above the value expected from the gas accre- 
tion rate, scaled by the dust-to-gas ratio. These values allow 
a dust trap such as a vortex or zonal flow to achieve a lo- 
cally critical over-density on the order of only a few tens of 
orbits and could trigger the streaming instability and gravo- 
turbulent formation of planetesimals, as discussed in Klahr 
& Bodenheimer (2006) and Johansen et al. (2006). 
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Appendix A: Simplified dust transport equation 

Assuming the gas to be the dynamically dominant medium, the 
radial transport of each trace species can be described by an 
advection-diffusion equation. 



dY.i 1 d 
dt r dr 



r \ Si Uj - Dj Sg — 







(A.l) 



where S, is the surface density of the tracer, Sg the gas surface 
density and m, and D,- are the tracer velocity and diffusivity, re- 
spectively. 

In our case, the two trace species considered are the small 
and the large representative grain sizes, as described in Section 4. 
The fraction of the total mass for each species is defined in 
Eq. 22. Thus, the evolution of the total dust surface density can 
be written as 



52d 

dt 



dt dt 
1 3 

K^oMo + Si Ml)] 

r or 



(A.2) 



1 d 
r dr 



SgDo 



dr 



+ 2gZ)i 



dr 



Using Eq. 22, we can simplify this to an advection-diffusion 
equation for the total dust surface density 



ffid 1 d 
£7? r or 



5:dM*-D*Sg- 
dr 



0, 



where 



u-{Di- D2) 



dfm{r) 



D* = {Di - D2) ■ f„,{r) + D2, 



(A.3) 



(A.4) 



(A.5) 



and u is the mass weighted velocity given by Eq. 24. The dust 
diffusivities are (see Youdin & Lithwick 2007) 



Di 



l+Stf 



(A.6) 



For Stokes numbers smaller than one, which is always fulfilled 
in this paper, the diffusivities are basically equal to the gas dif- 
fusivity and thus D* ~ D[ ~ D2 and u* ~ u. This allows us to 
simulate the radial evolution of the total dust surface density by 
using just the gas diffusivity and the mass weighted velocity, as 
in Eq. 23. 

Appendix B: Algorithm of the two-population model 

In this section, we summarize the algorithm of the two- 
population model. 

1 . calculate the representative size for a fragmentation limited 
size distribution afiag given by Eq. 8. 

2. calculate the representative size for a drift limited size distri- 
bution adrift given by Eq. 18. 

3. calculate the representative size in the case of drift induced 
fragmentation 



fldf = Stdf ■ 



2Eg 



(B.l) 



where Stdf given by Eq. 20. 
4. the representative size of the large population is now given 
by the smallest of the size limits 



ai = min(fldrift,afrag,fldf) 



(B.2) 



5. the initial phase where particles grow from the smallest size 
do to a 1 is approximated by 



fli(f) 



min 



«!, flo • exp 



t-tp 



where 



(B.3) 



(B.4) 



d i^k 



6. the velocities uq and mi of the two populations are then given 
by 



"0/1 



1 +Sti 



0/1 



where 



StO/l + Sty^j 

dlnP 



Mdrift, 



(B.5) 



"drift 



2Vi, d\nr 

is the drift velocity, P the gas pressure at the mid-plane 



(B.6) 



■Pg 



V2n 



and St the Stokes number is defined by 

flo/iPs n 
2' 



St, 



■0/1 



(B.7) 



(B.8) 



7. the effective dust transport velocity is the mass averaged ve- 
locity of both populations 



(1 - /m) • Mo + /m • Ml, 



(B.9) 



where the mass fraction depends on the size limiting 
mechanism 



fm 



0.97 if Adrift is the smallest of a^ag, adrift, ^df 
0.75 otherwise 

(B.IO) 

8. Having calculated the effective transport velocity m, the evo- 
lution of the dust surface density can be simulated by numer- 
ically solving the advection-diffusion equation, Eq. 23. 
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